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Recent microwave conductivity measurements of detwinned, high-purity, slightly overdoped 
YBa2Cu30e.993 crystals reveal a linear temperature dependence and a near-Drude lineshape for 
temperatures between 1 and 20 K and frequencies ranging from 1 to 75 GHz. Prior theoretical work 
has shown that simple models of scattering by point defects (impurities) in d-wave superconductors 
are inconsistent with these results. It has therefore been suggested that scattering by extended de- 
fects such as twin boundary remnants, left over from the detwinning process, may also be important. 
We calculate the self-energy and microwave conductivity in the self-consistent Born approximation 
(including vertex corrections) for a d-wave superconductor in the presence of scattering from ex- 
tended linear defects. We find that in the experimentally relevant limit (Q, 1/r <T< Ao), the 
resulting microwave conductivity has a linear temperature dependence and a near-Drude frequency 
dependence that agrees well with experiment. 



PACS numbers: 74.72.Bk, 72.10.Fk, 61.72.Mm, 78.70.Gq 



I. INTRODUCTION 

Of the numerous phases of the high-T c cuprate phase 
diagram, the low temperature superconducting phase ap- 
pears to be the least mysterious. All indications point 
to an order parameter of d x i_yi symmetry and well 
fined quasiparticle excitations above the condensate 
An excellent probe of the low energy excitations of this d- 
wave superconductor is the measurement of low temper- 
ature microwave conductivity. Experimentally, the real 
part of the microwave conductivity, cr(fl,T), can be ex- 
tracted from microwave-measurements of complex surface 
impedance Z s (£l, T)n~U Theoretically, at low tempera- 
tures, it is determined by calculating the linear response 
of a G?-wave superconductor in the presence of scattering 
due to static disorder J3~LL3 (Low temperature is-defined 
as low enough such that all inelastic scatteringEj is neg- 
ligible.) The mystery is that, in this least mysterious of 
phases, the simplest theoretical models do not agree with 
experiment £3 

The microwave conductivity of detwinned, high-purity, 
slightly overdoped YBa2Cu3 06.993 was measured by Hos- 
seini et alB Data was taken at five frequencies between 
1 and 75 GHz for temperatures from 1 to 95 K. Below 
20 K, in the regime dominated by elastic scattering, they 
find that the microwave conductivity has an approxi- 
mately linear temperature dependence and a near-Drude 
frequency dependence. Fitting the data to a Drude form, 
they extract a spectral weight that is linear in temper- 
ature (in agreement with the measured superfluid den- 
sity) and an effective scattering rate, 1/r, that is approx- 
imately constant over the entire low temperature regime 
(lto20K). 

From the extracted value of the scattering rate (5.6 x 
10 10 s _1 ) and the frequency and temperature range of 
the experiments, it can be seen that these measurements 
correspond to the parameter regime ft, 1/t <C T <C A 
where Aq is the gap maximum. For this parameter 



regime, the theoretical picture is as follows. A d-wave 
superconductor is characterized by an order parameter 
that vanishes linearly at each of four gap nodes. For 
T -c Ao, transport is dominated by low energy quasipar- 
ticle excitations generated in the vicinity of the nodes. 
For 1/t<T, the quasiparticles are generated thermally, 
rather than by the presence of disorder .113 (Hence, in what 
follows, we refer to this as the thermal regime.) The 
most straightforward model for such a system is one in 
which thermally generated nodal quasiparticles are scat- 
tered due to the presence of point defects such as impu- 
rities. Performing an impurity average and using a Kubo 
formula approach within the self-consistent t-matrix ap- 
proximation, the self-energy and. electrical conductivity 
have previously been calculated.ErHlill 

Surprisingly, the predictions of such calculations do not 
agree with the experimental results. The trouble comes 
in attempting to reproduce the constant scattering rate 
extracted from experiment. For a d- wave superconductor 
in the thermal regime, the anisotropic Dirac excitation 
spectrum yields a quasiparticle density of states that is 
linear in energy. This strong energy dependence of the 
density of states is then directly reflected in the results 
of the impurity t-matrix calculations. Such calculations 
yield 1/r ~ T and er ~ const in the Born limit and 1/r ~ 
1/T and a ~ T 2 in the unitary limit. More generally, 
over the full range of scattering strengths, it has been 
shown that simple modelstj of point defcqt-|(impurity) 
scattering are inconsistent with experiment .113 

Despite this, we certainly expect the presence of impu- 
rities to make a significant contribution to the total scat- 
tering. However, the observed disagreement between the- 
ory and experiment suggests that an additional scatter- 
ing mechanism may also be important. In what follows, 
we propose that this additional mechanism is scattering 
from extended linear defects. Our motivation for doing 
so is twofold. First of all, we note that scattering from a 
series of parallel lines differs from point defect scattering 
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in an important way. While point defects scatter in all di- 
rections, line defects only scatter in the direction normal 
to their orientation. Hence, whereas point defect scatter- 
ing samples the full quasiparticle density of states, line 
defect scattering restricts final states to those for which 
there is no change in the parallel component of momen- 
tum. Therefore, the line defect scattering rate need not 
inherit the strong energy dependence of the density of 
states and is thereby capable of exhibiting the less en- 
ergy dependent form required to agree with experiment. 
Secondly, we note that upon growth, YBCO crystals are 
full of "lines", twin boundary lines. Although the twin 
boundaries are subsequently removed via a dctwinning 
procedure, remnants of their structure may be left be- 
hind in the form of extended linear defects. Since lines 
are prevalent in this material and line defects tend to ex- 
hibit features of the measured behavior, it is logical that 
we should examine the contribution made by line defect 
scattering. For simplicity, we consider herein a system 
where extended linear defects provide the sole scattering 
mechanism, leaving the analysis of a system with both 
line defects and point defects for future study. 
In Sec. H 

we discuss twinning in YBCO and the na- 
ture of the detwinning process. In light of the observed 
twinning structure and guided by the notion that the de- 
twinning process leaves behind r emn ants in the form of 
line defects, we develop, in Sec. 
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a model for quasi- 
particle scattering from extended linear defects in a d- 
wave superconductor. In Sees. IV and |v|, we calculate 
the resulting self-energy and microwave conductivity. For 
parameter values within the thermal regime, we obtain 
analytical results which are compared to experiment in 
Sec. VIA. Numerical results, including deviations from 



our analytical expressions and valid beyond the thermal 
regime, are presented in Sec. VI B Conclusions are dis- 



cussed in Sec. VII where we provide physical motivation 
for our calculated results and discuss the implications of 
our findings. 



II. TWIN BOUNDARIES IN YBCOlf] 

Twinning occurs naturally during the growth of oxygen 
deficient YBa2Cu307_a. At full oxygenation (5 — 0), the 
crystal structure is orthorhombic at room temperature 
with lattice parameters a = 3.227 A and b = 3.8872 A, 
the latter being the orientation parallel to the CuO chain 
layer .EU For nonzero S, oxygen vacancies are present in 
the CuO layer and the remaining oxygen can disorder 
onto sites perpendicular to the CuO-chains, thereby re- 
ducing the orthorhombic distortion.E2l However, occupa- 
tion of these off-chain sites is inhibited by repulsive inter- 
actions with oxygen located at on-chain sites. The forma- 
tion of twin boundaries, which are abrupt boundaries be- 
tween crystal domains with CuO chain orientations that 
differ by 90°, are the energetic compromise between de- 
creasing the orthorhombicity and maintaining the CuO 




FIG. 1. Photograph of as-grown (twinned) and detwinned 
single crystals of YBCO as viewed through an optical mi- 
croscope under polarized light. In the as-grown sample, two 
different CuO chain orientations, visible as a difference in 
grayscale, are separated by twin boundaries oriented at ±45° 
to the horizontal. (Arrows highlight the two CuO chain ori- 
entations.) Darkened areas are regions of highly concentrated 
twin boundaries. In the detwinned sample, twin boundaries 
are no longer visible and the a axis is parallel to the horizontal 
throughout the crystal. (Photograph courtesy of R. Harris.tj) 



chain ordering. 

The overdoped YBa2Cu306.gg3 crystal studied by Hos 
seini et at was prepared by annealing an as-grown cry; 
tal at 350° C for 50 days in flowing high purity oxygen 
Detwinning was accomplished by incrementally apply- 
ing uniaxial stress in the ab plane with the sample tem- 
perature fixed at 250°C in flowing high purity oxygen. 
Progress was monitored by viewing the sample through 
a microscope objective with a polarizer oriented at 45° 
with respect to the crystalline axes. Domains with dif- 
fering CuO chain orientations were then visible and the 
elimination of twin boundaries could be verified. The 
complete detwinning of a 1 mm 2 crystal can typically be 
accomplished within one day. 

An example of an as-grown twinned YBa2Cu307_^ 
sample and a detwinned sample of YBa2Cri306.gg3 are 
shown in Fig. [l| This picture was obtained by shining 
polarized white light onto the crystals at 45° with re- 
spect to the crystalline axes and then observing through 
a polarizer oriented parallel to the incident light. In the 
as-grown sample, the two different CuO chain orienta- 
tions are visible as a difference in grayscale. The lines 
oriented at ±45° which separate these regions are the 
twin boundaries. Darkened areas denote regions of highly 
concentrated twin boundaries. In the detwinned sample, 
there is a single orientation for all the CuO chains. Here, 



2 



the a axis is parallel to the horizontal since this is the 
direction along which the uniaxial stress was applied. 
Note the absence of any visible twin boundaries in the 
detwinned crystal. 

III. EXTENDED LINEAR DEFECT 
SCATTERING 

Since line defect scattering has the appealing property 
that it does not sample the full quasiparticle density of 
states, we consider a model in which d-wave quasiparti- 
cles scatter from extended linear defects. In particular, 
we imagine that the line defects are arranged in a man- 
ner reflective of the twin boundary structure exhibited in 
YBCO prior to detwinning. That is, we consider domains 
of line defects within which the lines are parallel, oriented 
at either +45° or —45° to the crystal axes, and separated 
by distances on the order of microns. Our picture is that 
in the process of detwinning, neighboring twin bound- 
aries are made to annihilate in order to eliminate regions 
in which the CuO chain orientation is disfavored by the 
applied uniaxial stress. Although this procedure is effec- 
tive in removing the twin boundaries, it can leave behind 
defects along the lines at which the twin boundaries an- 
nihilate. While there may be multiple ways in which this 
can come to be, one example of a detwinning scenario 
through which line defects could be left behind is the 
following. 

The twin boundaries in as-grown YBCO are ±45° lines 
which separate regions with horizontal CuO chains from 
regions with vertical CuO chains. Geometry dictates that 
the distance between oxygen sites on opposite sides of 
such a line is smaller than the normal oxygen-oxygen 
separation in the bulk. As a result, it is energetically ad- 
vantageous for oxygen vacancies, which must be present 
in the CuO layer of doped YBCO, to be concentrated 
along twin boundary lines. Upon detwinning, a uniax- 
ial stress is applied to favor one CuO chain orientation 
and force neighboring twin boundaries to approach each 
other and annihilate. As the twin boundaries move, the 
vacancies move with them until neighboring boundaries 
annihilate, leaving behind lines of oxygen vacancies. The 
resulting potential, felt by quasiparticles in the Cu02 
plane, is that of a collection of extended linear defects 
arranged in a pattern that reflects the original twinning 
structure. Note that this is just one example of a process 
that could yield line defects. As we shall see, however, 
our results will suggest that something similar to this is 
going on. 

Consider a single domain in which all line defects are 
aligned at —45° to the horizontal axis. This situation 
is depicted schematically in Fig. |[ We denote the di- 
rections parallel and perpendicular to the defect lines as 
f || and fj_ respectively. Now consider the effect of scat- 
tering off of such extended linear defects in 2d. Since 
all the lines are normal to the rj_ direction, they act as 




FIG. 2. Schematic depiction of a single domain of extended 
linear defects. The line defects are aligned at —45° to the 
horizontal and spaced randomly in the fx direction with a Id 
density nr- 

potential sources, V(r±), that depend only on the r± 
coordinate. Fourier transforming, it is clear that such 
scattering events conserve the parallel component of mo- 
mentum, py, and are therefore one-dimensional in nature. 
Furthermore, since the distribution of line defects is ran- 
dom in the rj_ direction, this problem corresponds to an 
effective one-dimensional "impurity" problem (along f j_ ) 
in which p|| is conserved but not relevant to the scat- 
tering. Performing a Id disorder average, properties of 
the resulting system can be calculated in terms of the Id 
line defect density, nx, and the Id momentum-space line 
defect scattering potential, V(p± — k±). 

Now consider the effect of this sort of one-dimensional 
scattering on the physics of a d-wave superconductor. We 
model a generic d-wave superconductor as a system with 
the Brillouin zone of a 2d square lattice and an order 
parameter of d x 2_ y 2 symmetry that vanishes at each of 
four nodal points on the Fermi surface. In the vicinity 
of each of these nodal points, the electronic dispersion, 
ek , varies linearly across the Fermi surface and the order 
parameter, varies linearly along the Fermi surface. 
As a result, near each of the gap nodes, the Bogoliubov 
quasiparticle excitation spectrum takes the anisotropic 

Dirac form E\. = \J e\ + Al = \Jv 2 fk\ + v\k\ where the 

degree of anisotropy is expressed by the ratio of the Fermi 
velocity, Vf, to the gap velocity (slope), «2, and the mo- 
menta, k\ and k2, are defined locally at each node such 
that in all cases, ej; = v/fci and Aj, = t^fe. Hence, 
excitations are free at the nodal points and at temper- 
atures much less than the gap maximum, quasiparticles 
are thermally generated within small regions about the 
gap nodes. The situation is depicted in Fig. |^ where the 
elliptical regions (labeled 1 through 4) denote four pock- 
ets of thermally generated quasiparticles. In the case of 
point defect (impurity) scattering, quasiparticles could 
be scattered either within the nodal pocket from which 
they originated (intra-node) or from one node to another 
(inter- node) .E3 This is so because point defect scattering 
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FIG. 3. Schematic depiction of the Brillouin zone of a 
d-wave superconductor in the presence of line defect scatter- 
ing. The elliptical regions, labeled 1 through 4, denote (at a 
very exaggerated scale) four pockets of thermally generated 
nodal quasiparticles. In addition to the global momentum di- 
rections, p|| and kx, defined parallel and perpendicular to the 
defect lines, we have also defined local momentum directions, 
ki and k2, about each of the four nodes. The dashed lines 
denote scattering paths that conserve the parallel component 
of momentum. 

events do not conserve either component of momentum. 
However, due to the Id nature of extended linear defect 
scattering, scattering events conserve the parallel com- 
ponent of momentum, pu . Therefore, in the presence of 
line defect scattering, quasiparticles can only be scattered 
along lines parallel to the k± axis in momentum space. 
Four such (dashed) lines are shown in Fig. From the 
figure it is clear that the condition of p» conservation 
means that we have two different kinds of nodes. For 
the odd nodes (1 and 3), quasiparticles can be scattered 
either intra-node (1 — > 1 or 3 — > 3) or to their opposite- 
node (1 — > 3 or 3 — > 1). Forbidden by pii conservation is 
the adjacent-node scattering (i.e. 1 — ► 2) that would be 
allowed for the point defect case. But for the even nodes 
(2 and 4) , quasiparticles can only be scattered intra-node 
(2 — > 2 or 4 — > 4). Here both adjacent-node and opposite- 
node scattering are forbidden by pn conservation. Due to 
the clear difference between odd nodes and even nodes 
in the presence of line defect scattering, it will be impor- 
tant in the calculations that follow to always treat the 
odd node and even node cases separately. Furthermore, 
we should note that the designation of odd or even to a 
particular node is a direct result of our choice to consider 
a domain in which the line defects arc aligned at —45° 
to the horizontal. To consider the other type of domain, 
where the line defects are aligned at +45° to the horizon- 
tal, we need only swap odd for even in all designations. 
With this model of thermally excited pockets of quasipar- 
ticles scattered by parallel randomly-spaced defect lines 
in mind, we proceed to calculate the resulting self-energy 
and microwave conductivity. 



IV. SELF-ENERGY CALCULATION 

For a superconductor at finite temperature, our cal- 
culations will employ Matsubara Green's functions ex- 
pressed in the 2x2 matrix Nambu formalism. The bare 
Green's function takes the form 

where the tilde denotes a Nambu-space matrix and iuj — 
i(2n+l)irkBT is a fermionic Matsubara frequency. In the 
presence of scattering, the Green's function is dressed via 
Dyson's equation 

£?(k, ilu)- 1 = g {k, iuj)- 1 - E(k, iuj) (4.2) 

where 



S(k, ilu) — S(fcj_, fen , iu>) = 




is the Matsubara self-energy. 

In the Born approximation, this self-energy matrix can 
be calculated by evaluating the diagram in Fig. | to ob- 
tain 

Z{p±,p\\,iuj) = n T y~] \V(px - k±)\ 2 f 3 G(kx,p\\,iu)f3 

(4.4) 

where V(p± — fej_) is the line defect scattering potential, 
tit is the line defect density, and the fj are Pauli matrices 
in Nambu-space. Due to the Id nature of the line defect 
scattering, pu is conserved and we have integrated only 
over fcj_ . Since quasiparticles reside only in the vicinity of 
the four gap nodes, each momentum can be expressed by 
a node index, j, and a local momentum, (ki, k 2 ), defined 
about node j. Furthermore, it is useful to scale out the 
anisotropy of the excitation spectrum by defining local 
scaled momenta, k[ = Vfki and k' 2 = ^2^2, such that 
at each node ek = k[, = k' 2 , and Ek = y kf + k' 2 2 . 
For convenience, we drop the primes in all that follows 
and take all locally defined momenta to be scaled mo- 
menta. Making use of these locally defined variables, we 
can replace 

E - E /°° ^ («) 



Vfpj-kj)/ \V(k - p ) 




FIG. 4. Self-energy within Born approximation. 
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where the sum over node index j is restricted to include 
only nodes that are crossed by the k± integration path, k 
is the local scaled momentum component parallel to k± 
(k\ or k 2 ), and v is the corresponding velocity (vf or v 2 ). 
Similarly, we can replace the initial momentum, (p_L,pii), 
with an initial node index and an initial local scaled mo- 
mentum, (jpiiPi). As is clear from the dashed integration 
paths shown in Fig. |[ the details of these replacements 
depend on the type of node about which the quasiparti- 
cles reside prior to scattering. For quasiparticles initially 
near odd nodes (1 or 3), the node index sum yields both 
an intra-node scattering term and an opposite-node scat- 
tering term, in both of which k\ is the local integration 
variable and p 2 is the local conserved variable. In con- 
trast, for quasiparticles initially near even nodes (2 or 4), 
the node index sum yields only an intra-node scattering 
term, for which k 2 is the local integration variable and 
Pi is the local conserved variable. It is therefore clear 
that onc e we change to local coordinates, we obtain from 
Eq. (4.4) two self-energy expressions, one for odd nodes 
and one for even nodes. For odd nodes: 

n T = f dki \ rr2flol 



E°(p2, iw) = 77—^3 / — V 1 'g°(k 1 ,p 2 ,iu) 
2Vf J 7T 



+W(* l ,-P2,«o) r 3 (4.6) 



g j (qi,q 2 ,iuj) 



(im - + (ft + Sj)f 3 + (q 2 + g|)n 

(<?2 + S J 2 ) 2 
(4.9) 



(iw - SJ) 2 - ( qi + E{) 2 



where (qi,q 2 ) is the scaled momentum about a particu- 
lar node and the self-energy components are functions of 
q 2 for odd nodes and q\ for even nodes. This expression 
for the Green's functions in term s of the self-energies to- 
gether with Eqs. (4.6) and ( |4.7D , giving the self-energies 
in terms of the Green's functions, comprise matrix equa- 
tions which can be solved self-consistently for the odd- 
node and even-node self-energies. 

Let us consider the odd case first. Shifting the inte- 
gration variable k\ by E°, we see that the f 3 term in- 
tegrates to zero and therefore EJ = 0. Then continuing 
iuj — ► U3 + iS and making the ansatz (which will soon be 
shown valid) that E° is an even function of p 2 while E£ 
is an odd function of p 2 , we obtain 



E° 



£2 



dk\ 



a°{io - E°) 



7T (w - £°) 2 - (pa + E°) 2 - k\ 



dki 



-a?(Pa 



7T (ijj ~ E°) 2 - (p 2 + S^) 2 — kf 



(4.10a) 



(4.10b) 



whereas for even nodes: 



E (P1,IUJ) = —T 3 

2v 2 



dk 2 



VtQ e (j>iMM f 3 (4.7) 



where the superscripts 'o' and 'e' denote odd and even 
respectively. Here we have parameterized the scattering 
potential using the notation defined in Ref. whereby 
V\ defines the potential for an intra-node scattering event 
and V3 defines the potential for an opposite-node scatter- 
ing event. Note that for both odd nodes and even nodes, 
the self-energy is a function of frequency as well as the 
conserved component of the local scaled momentum, p 2 
for the odd case and p\ for the even case, but never the 
momentum component along k±. This is a tremendous 
simplification since it means that the sel f-en ergy func- 
tions built into the right-hand side of Eqs. (4.6) and (4.7) 
will not be functions of the integration variables and can 
always be treated as constants with respect to the inte- 
grals. It is this fact that makes these equations tractable. 

Since the bare Green's function, Eq. (4.1), has no r 2 
com ponent, it is clear from the form of Eqs. (4.2) and 
(4.4) that it is self-consistent for the self-energy matrix 
to lack a f 2 component as well. It is therefore convenient 
to write the odd-node and even-node self-energy matrices 
in the form 



(4.8) 



where j = {o, e}. It then follows from Eq. (4.2) that the 
corresponding dressed Green's functions can be expressed 



as 



n T {V? 



Vf) 



n T {V? 



Vi) 



2 V/ 



2v f 



(4.10c) 



where S° = Y^°(p 2 ,u>) and S 2 = T,2(p 2 ,uj) are now re- 
tarded functions. 

For even nodes, we follow an analogous procedure 
and obtain results of the same form. This time, shift- 
ing the integration variable k 2 by S|, it is the fx term 
that integrates to zero such that T, 2 = 0. Continuing 
iuj — > lu + iS to obtain retarded functions S e = S e {p\ , uj) 
and Sf = Sf (pi , u>) , we find 



E e = 



dk 9 



a e (uj - S e ) 



7T (w - E e ) 2 - (pi + Ef 



e\1 



k 2 



dko 



af(Pi 



7T (w - E e ) 2 - (pi + Ef 



e ri T V x 2 
a = a, 



e\ 2 



k\ 



2v 2 



(4.11a) 



(4.11b) 



(4.11c) 



where we have defined a\ in analogy with the odd node 
case. 

In light of the similarities between the odd and even 
results, we can treat both cases on the same footing by 
defining for j = {odd, even}: k — {ki,k 2 }, p = {p2,Pi}, 
E = {S°,E e }, Ei = {E 2 , E^}, a = {a°,a e }, qj = 
{a°,al}, and i] = {— 1, +1}. Hence, keeping in mind 
that our new variables have different meanings for odd 
and even nodes, we can write 
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e = 



dk 



7T (w-S) 2 -(p+S 1 ) 2 -fc 2 



(4.12a) 



Si 



r/ai(p + Si) 



7T ( W -S)2_(p +Sl )S 



fc 2 



Tj- 



a.\ p + Si , 



(4.12b) 



where the second eq uality in Eq. (4.12b) w as ach ieved by 
making use of Eq. ( [4.12a]) - Solving Eq. ( |4.12b|) f or Ei 
in terms of E, plugging the result into Eq. (4.12a), and 
evaluating the integral we find 



(w-7e) 



1/2 



Ei 



(7- 

LU — 7E 



(4.13) 



where we have defined 7 = 1 + rjai/a . This equation will 
be solved numerically in Sec. VI B to obtain an exact 
result for the self-energy. 

Fortunately, in the experimentally relevant "thermal 
regime", ft, 1/r -C T <C Aq, we can extract an approxi- 
mate analytic result. As will be shown in the microwave 
conductivity calculation discussed in the following sec- 
tion, the energies of interest are |w| on the order of T 
and what we have called a is on the order of 1/r. Thus, 
in the thermal regime, we are justified in taking the limit 
|w| ^> a such that we can set uo — 7E m u. As a result, 
defining real and imaginary parts via E = A — iT and 
Ei = Ai — iT\, we find the approximate expressions 



^■-W-W) A 



1- * 



—a sgn(w) 



6Qp\-\w\) 



(S) 2 -l 



(4.14a) 



Ti w ?yai 



p 0(M - \p\) 



Ai 



-■qai- 



P 0(\P\ ~ M) 



1- * 



^ -1 



(4.14b) 



where the sharp cutoffs of the theta-functions are clearly 
artifacts of our approximation that are smoothed away 
in the exact solution. In the following section, we will 
make use of these thermal-regime self-energy functions to 
obtain an analytic expression for the microwave conduc- 
tivity in the experimentally relevant parameter regime. 



V. MICROWAVE CONDUCTIVITY 
CALCULATION 

The electrical conductivity tensor can be calculated by 
means of the Kubo formula 



a(ft,T) = - 



ImII ret (ft) 
ft 



(5.1) 




(b) > 

A 




FIG. 5. Polarization bubble and dressed vertex within 
Born approximation. 



where n rei (ft) = II(ift -> ft + i8) and Il(ift) is the Mat- 
subara polarization function (or current-current correla- 
tion function). Evaluating the diagram in Fig. |5|(a) we 
find 

%uj p 

x Tr \g(p, iuj)G{p, iuj + ift)f (p, iu, ift) (5.2) 

where Vf p points in the direction of the Fermi velocity 
and r is the dressed vertex function. Including vertex 
corrections within the Born approximation, T is calcu- 
lated by evaluating the ladder diagrams in Fig. ||(b). Do- 
ing so, we obtain 



r( P j 



iu, id) = v/pi + ^2,tit \V(p± - k±)\' 



xf 3 G(k±,p\\,iuj + id)t(k±,p\\,iu,id)g(k±,p\\,iu)f 3 

(5.3) 

where tit is the line defect density and V(p± —k±) is the 
line defect scattering potential. Once again, due to the 
Id nature of our scattering, the momentum component 
parallel to the defect lines, p\\, is conserved and we only 
integrate over the perpendicular component, k±. Replac- 
ing momentum integrals by sums over node index and in- 
tegrals over local scaled momenta, the polarization tensor 
takes the form 



j=l P iu, J 



d 2 p 



(2-K) 2 VfV2 



xTr 



G j (p, iu))g j (p, iu + id) (i + A J (p, iu, id) 



(5.4) 



where we have defined a vertex correction function, K° , 
such that 



f 3 (p, iu,id) = v j f (i + A J (p, iu, id) 



(5.5) 



Note that A J depends only on the conserved component 
of the scaled local momentum (i.e. A° = h°(p2,iu,id) 



G 



while A e = A e (pi, iui, in)). Since quasiparticles in odd 
nodes (1 or 3) can be scattered either intra- node or to 
the opposite-node while quasiparticles in even nodes (2 
or 4) can only be scattered intra-node, we must treat the 
odd and even cases separately as we proceed to calculate 
the vertex correction. 



For odd nodes, Eqs. (5.3) and ( |5.5| ) yield 

2v f J 7T 



(i + A°( P2 ))g (k 1 ,p2,iuj)f3 

[/? r die., 

f 3 G°{kx, —p 2 , + in) 



nrVj f dki . fio 
2v f 



x (t +A o (-p 2 ))0 o (fc 1 ,-p 2 ,iw)f 3 



(5.6) 



where we have suppressed the frequency dependence of 
A° for simplicity. Taking the Nambu-space trace, noting 
that f 3 G°(p 2 )f 3 = G°(—p 2 ), and using the cyclic proper- 
ties of the trace, we find 



Tr 



A°(p 2 ) 



= Tr 



2v f 

+VfA°(jp 2 ) - V 3 2 f 3 A°(-p 2 )f 3 



(5.7) 



where we have defined 

i°{P2) = 

= r{ P2 )t + r 1 {p 2 )T 1 



dki ~ 

G°(kx,p 2 , iuj)g°(ki,p 2 ,iuj + iti) 

ir 



(5.8) 



with I°(p 2 ) an even function and I°{p 2 ) an odd function. 
Note that the fact that I°{p 2 ) can be written in this 
form is merely assumed at this point but will be demon- 
strated later. Manipulating further and making the addi- 
tional assumption (which will soon be shown valid) that 
Tr[A°(p 2 )] is an even function and Tr[fiA°(p 2 )] is an odd 
function, we can write 



Tr 



A° 



Tr 



A" 



a\r x Tr 



fiA° 



(5.9) 



where a° = n T (V 2 - V 3 2 )/2v f , as defined in Eq. ( |4.10c[ ), 
and we have now suppressed the p 2 dependences. Going 
back to Eq. ( |5.6| ) but this time multiplying by fx before 
taking the trace, we can proceed along analogous lines to 
show that 



Tr 



fxA° 



= a° 1 J°Tx 



fxA° 



a\Jl 2 + Tr 



A" 



(5.10) 



where we have now defined 

J°(P2) = 

= J?(p 2 )l + J°{p 2 )fx 



dki ~ 

§°(kx,p 2 , iu)fxG°(kx,p 2 ,iw + in) 

IT 



(5.11) 



with J°{p 2 ) a n ev en fu nction and J°(p 2 ) an odd function. 
Solving Eqs. (5.E) and ( |5.10 ) simultaneously, we find that 



Tr 



A" 



2a<j (iQ + g&L) 



(5.12) 



and 



Tr 



fxA° 



1-q°J° 



(5.13) 



from which it is clear, given the defined parity of 1°, 
Ii, J°, and J°, that as a function of p 2 , Tr[A°] is odd 
and Tr[fiA°] is even, in agreement our prior assumptions. 
Finally, using Eqs. (^) and ( |5.9| ) and the fact that the 
Pauli matrices are traceless, we can write 



Tr 



A°( P2 ) 



a?Tr 



I°(p 2 ) (l + A°(p 2 ) 



(5.14) 



which will prove a useful expression in what follows. 

For even nodes, the vertex correction function takes 
the form 



A e (Pi) 



n T V? 
2v 2 



dk 



2 T 3 G e {pi,k 2 ,iuj + in) 



x (l+A e {px)) G e (px,k 2 Mh- 



(5.15) 



Note that this even-node expression is of the s ame form as 
the corresponding odd-node expression, Eq. (5.6), but is 
somewhat simpler since we need only consider intra-node 
scattering. Taking the trace and cyclically permuting 
within the trace, we find 



Tr 



A e ( Pl ) =alTr I e ( Pl ) 1 + A e ( Pl ) 



(5.16) 



where a\ — n^V 2 /2vf, as defined in Eq. (4.11c), and we 
have defined 



I e (Pi) 



^-G e (p 1: k 2 , iu)G e {pi, k 2 ,iuj + in) 

7T 



I e (px)t + It(px)T 3 . 



(5.17) 



Making use of the second equality in Eq. ( |5.17| ) and fur- 
ther defining 



J e (pi 



dk 2 ~ 



G e (px,k 2 ,iuj)T 3 g e (p 1 ,k 2 , ilu + in) 



= Jt(px)l + J e (pi)f 3 (5.18) 
we obtain two coupled equations for Tr[A e ] and Tr^A 6 ], 



Tr 



A e 



all e 2 + Tr 



A' 



a\If Tr 



(5.19) 
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Tr 



r 3 A c 



r 3 A £ 



t\Jl ( 



2 + Tr 



A' 



Solving simultaneously yields 



Tr 



A' 



1(1 



2a? I 



1-ctf J' 



1 



(5.20) 



(5.21) 



which has precisely the same form as in the odd case. 

Due to this similarity in form, we can unite the odd 
and even cases by using the superscript j — {o, e} and 
defining a generalized integration variable k — {ki,k 2 } 
and a generalized conserved variable p — {p 2l pi}. Doing 
so, we can define 

1 



2a{ 



K 3 (p, iuj, ity = — Tr A 3 {p, ito, in) 

p + AllA. 



a J V J 3 
1—cA Ji 



1 - a{ [P 

and note via Eqs. ( |5.14 ) and ( 5.16 ) that 



(5.22) 



Tr 



V (p, iuj, ift) ( 1 + A 3 (p, iu, in) 



2K j (p,iu,in). 

(5.23) 



Then going bac k to our expression for the polarization 
tensor, Eq. (|5.4|), and noting that 



v / v / 

« 2 - 2 
V / V / 



V / V / 



Tl 
Tl 



(5.24) 



we can write the tensor as the sum of an odd-node term 
and an even-node term 



IT(iQ) = n°(^) 



II 



n e {in) 



1 - Tl 



(5.25) 



Matsubara Green's function. A procedure for evaluating 
Matsubara sums of functions of this type has been devel- 
oped in Appendix B of Ref. Applying this procedure 
to the case at hand, we find that the imaginary part of 
the retarded polarization function takes the form 



2ir v 2 
xRe 



dp 



doj 
2^ 



(np(uj +0) — np(oj)) 



K B (p,u,n)-K A (p,u,n) (5.27) 



where rip is the Fermi function and 

K 3 A (p, oj, ft) = lim K 3 (p, oj + iS, itt) 
K j B (p,uj,n) = lim K j (p,(J-iS,iQ) (5.28) 

are what we will refer to as the A-form and B-form of 
K 3 (p, iu, ifl). Plugging into the Kubo formula, Eq. ([Tl|), 
this yields the conductivity tensor 

ct(0, T) = a°(n, T) \l + ri] + cr e (Sl, T) \t - ri 

(5.29) 

where 

/ riF (w) — ni?(w + n) * 



*>{n,T)- 



e 2 v f 



, duo 

4?H V 2 J_oo V 



x / dp Re 



K 3 B (p,w,Sl)-K 3 A (p,u,Qi) . (5.30) 



The next step is to express K 3 A and K B in terms 
of the self-energy calculated in Sec. To do so we 

must first evaluate our expressions for the Matsubara- 
formalism functions, P, If, J 3 , and J{, and then an- 
alytically continue to obtain the corresponding A-forms 
and B-forms. Dressed via the self-energy functions, the 
Matsubara Green's functions, evaluated at frequencies ito 
and iuj + in, take the form 



where, for j = {o, e}, 



—Q 3 (k,p,iuj) 

7T 



e 2 v f 



xQ 3 (p,k,iuj + iti) (i + A j (p,iw,iQ,) 

K j (p,iu,iCl). (5.26) 



2nv 2 J 



At this point, we would like to evaluate the Matsubara 
sum of K 3 (p,iuj,ifl), analytically continue the Matsub- 
ara frequencies, and obtain retarded functions. To do so 
we should note that the internal and external frequencies, 
iuj and in, enter K 3 (p,iuj,in) only through functions 
of "Matsubara-couplets" of the form A(iu>)B(iui + in), 
where both A and B have the analytic structure of a 



g°(k lP ,tLo) 



9f 



y (k,p, iuj + iil) = 



for odd nodes and 



G e (p,k,tu;) 



fr~k 2 -g° 2 



f e*_ g e2_ k2 



o2 



(5.31a) 



(5.31b) 



(5.31c) 



g e (p, k, iuj + in) = —2 ^ 7J- (5.31d) 

J 2 ~ 9 2 ~ fc 

for even nodes, where we have defined 
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fl=iu- E J (p, iw) 
/| = ioj + iO, — TP (p, iu + id) 
g{ = p + T,{(p,iw) 

Plugg in g the dress ed G reen's functions into Eqs. 
(|5 11[), (p.l7|), and (|5.18|), we find that 



(5.32a) 
(5.32b) 
(5.32c) 
(5.32d) 



/1/2 + 9x9% + k 2 



tt (A 2 



^ 2 )(/| 



<?1 



k 2 ) 

(5.33a) 



1\ (p, iu/, ifl) — r/Ji (p, iu, zf2) 

rffc fm + fi9\ 



* (fl-9!-k 2 )(fi-g 2 2 -k 2 ) 

(5.33b) 



J(p, iu, ifl) = i] 



dk 



.fih+9192 - k 2 



n (f 2 g\ k 2 )(f 2 g 2 k 2 ) 

(5.33c) 



where 77 = {— 1, 1} for j = {o, e} and we have suppressed 
the node-index superscripts for simplicity. Expanding 
the integrands via partial fraction decomposition, noting 
that the result ing in tegrals are precisely of the form that 
appear in Eq. ( 4.12 ), and using this to relate the integrals 
to our self-energy functions, the above can be e valua ted. 
Then analytically continuing as prescribed in (5.28), we 
obtain both the A-form functions 

_ l (/ 2 + A)(E - E+) + v £(g 2 - giKSi + E 1+ ) 

A ~ * (fi-f?)-(g 2 2-9i) 

(5.34a) 



l (g 2 E-g 1 S] + )+^(/ 2 S] 1 -/igi + ) 
Ha = 772 7T\ r~2 2^ (5.34b 

a {fi - ft) - (gi - at) 



J A = 



r, (/ a - /i)(E + S+) + V ^(g 2 + g 1 )(E 1 - E 1+ ) 



(fl-m-(9 2 2-9l) 



(5.34c) 



and the B-form functions 



_ 1 (h + A*)(E* - S+) + V^(92 - 9tm + Si+) 

(fi-ff)-(9 2 2 -gf) 



Q 



(5.35a) 



i(a a s*-5 1 *s+) + » 7 ^(/ a E;-/ 1 *E 1+ ) 

/lB ~ a (ft ~ I?) - (9t ~ 9f) (5 - 35b) 



(/ 2 -/ 1 *)(E* + E 4 



•^)(E*-E 1+ ) 



(/I - /1 



*2\ 



(9 2 2 



9f) 



(5.35c) 



where all self-energies are now retarded and E = E(p, w), 
E + = E(p,w + fi), Ei = Ei(p, w), Ei+ = Ei(p,w + fi), 
/1 = u - E, / 2 = u + - E + , g x = p + Ex, and g 2 = p + 
Ei + . Then, in terms of the A-form and B-form functions 
above, 



K, 



l(A,B) 



1-OLxJf. 



(A,B) 



Oil If 



(A,B) 



1 a ' J i%,B) 
1 — CIl J(A,B) 



(5.36) 



Thus , g iven t he se lf- energ y func tions , E and Ei, Eqs. 
( |5~29l ), ( ^30| ), ^3§, (|5T355 , and constitute a well- 

defined pres cripti on for calculating the conductivity ten- 
sor. In Sec. VI B, these equations will be evaluated nu- 



merically, in conjunction with a numerical solution for 
the self-energy, in order to obtain an exact result for the 
conductivity. 

However, in the thermal regime, Q, 1/r <C T <C Ao, 
we are justified in taking the limit Q, a,ot\ <C \u\, which 
simplifies our results significantly and allows for an an- 
alytic solution. Taking this limit and using the thermal 
regime self-energy calculated in Sec. |v|, we find that the 
B-form results become 



u 

Ib ~ —I\B 
P 



where 



' 1 P ' Ju 



6(\u\ 



1 + ~iw 
-\P\) 



1 



(5.37) 



(5.38) 



as expressed in Eq. ( 4.14 ). Remarkably, these functions 
conspire to give Kb the surprisingly simple form 



Kr 



2r tr 



1 — iflTtrY 



(5.39) 



where we have suggestively defined an effective transport 
scattering rate 



1 

2r tr 



OL\. 



(5.40) 



Calculating the A-form functions to the same order as in 
the the B-form case, we find that 



I A ~ I\A ~ J A ~ 



K A ^0. 



(5.41) 



Therefore, taking the real part of Kb — Ka and integrat- 
ing over scaled momentum we see that 



/oo 
dpRc[K B -K A } =4/(fiT tr )- 
-00 



Ttr 



M (5.42) 
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where we have defined the function 



/(*) 



arctanh 




(5.43) 



which goes as 1 — x 2 /3 for small argu ment and as ln(2x) 
for large argument. Plugging into Eq. (5.30), noting that 



dui \w\ 



(5.44) 



and writing the node-index superscripts explicitly, yields 
the node-j microwave conductivity 



^(fi,T) = 4- 21n ( 2 )/(^' ; " 



It is important to note that for odd nodes 



k B T. (5.45) 



1 _ n T {V 2 + V 2 ) n T {V 2 - V 2 ) _ n T V 2 



2t& 



2v f 



2v f 



while for even nodes 



1 n T V 2 n T V 2 



2tI 



2v 2 



2v 2 



Hence, for non-zero frequencies 



1 







<r e (f2,T) = 



(5.46) 



(5.47) 



(5.48) 



so the even nodes make no contribution to the microwave 
conductivity. From a physical point of view, these re- 
sults are quite logical. Since electrical current is propor- 
tional to Fermi velocity, it is directed outward along the 
Brillouin zone diagonals at each node (see Refs. p4| , p3"| ). 
Therefore, quasiparticles in the vicinity of odd nodes 
carry current normal to the defect lines which is de- 
graded via line defect scattering. However, quasiparticles 
about even nodes carry current parallel to the defect lines 
which is unaffected by their presence. Since a nonzero mi- 
crowave conductivity requires a nonzero scattering rate, 
only the odd nodes can contribute. Hence 



a(Cl,T) = <r°(n,T) 



a 



(5.49) 



which clearly contains off-diagonal components and is 
therefore spatially anisotropic. This is the case because 
we have been considering only a single domain in which 
all defect lines are aligned at —45° to the horizontal. 
Since such a system certainly has a preferred direction, 
an anisotropic conductivity makes sense. However, in re- 
ality, we expect a detwinned crystal to have both —45° 
and +45° aligned domains. As discussed in Sec. Ill, our 
results can be adapted to the +45° case by swapping all 
odd designations for even and vice versa . Doi ng so simply 
changes the sign of the *T\ term in Eq. ( |5. 49 ). Therefore, 



averaging over domains with ±45° defect lines cancels 
the off-diagonal components and yields an isotropic con- 
ductivity tensor 



e 2 v f 



i?(n,T) = —-L2]n{2)f(Qr t ° r ) r 



: k B Tl. (5.50) 



Hence, within the thermal regime, Born scattering from 
extended linear defects in a (i-wave superconductor yields 
both a linear temperature dependence and a near-Drudc 
frequency dependence. 



VI. RESULTS 

A. Analytical Results 

Now let us compare our thermal regime expression for 
the microwave conductivity due to extended linear de- 
fect scattering, Eq. (5.50), with the measured microwave 
conductivity of detwinned YRa^CuaOg.ggs obtained via 
experiment by Hosseini et alB A fit of our calculated 
result to the temperature-dependent part of the mea- 
sured conductivity data is presented in Fig. Note 
that the close agreement seen in the figure was achieved 
with only two free parameters: the effective transport 
lifetime, r tr , and an overall scale factor. The fit yields 
a lifetime of 2.93 x 10~ n s which corresponds to a mean 
free path on the order of microns (quite reasonable for the 
high-purity sample in question) . Assuming an anisotropy 
ratio, Vf/v2 — 21, as measured via sub-.Kelvin thermal 
conductivity by Taillefer and co-workersEj, we obtain a 
scale factor of 0.6 by which our expression must be mul- 
tiplied to fit the data. This isjp-Lip^ with the expected 



Fermi-liquid correction factorc3c3c3, a 2 ^ ~ 0.4 — 0.5 
that has beep-^J^taiped from measurements of the super- 
fluid densityuO'D'El Thus, our thermal regime expression 
yields a quantitative fit to the temperature-dependent 
part of the measured data. 

However, there are several featur es o f the Hosseini 
data that are not captured by Eq. ( 5.50 ). First of all, 
for each of the frequencies at which data was taken, a 
temperature-independent shift was observed in addition 
to the temperature-dependent part plotted in Fig. ^. Yet 
in our thermal regime expression, cr(£l, T = 0) is zero. 
Furthermore, though predominantly linear with tempera- 
ture, the measured data deviates slightly, but noticeably, 
from linearity. In fact, Hosseini et al. note a gradual 
evolution from a concave-down (sub-linear) deviation at 
low frequencies to a concave-up (super-linear) deviation 
at high frequencies. Our thermal regime expression is 
strictly linear with temperature. 

From the measured data (see Fig. 4 of Ref. ||), it 
seems clear that the observed concave-down deviation 
at low frequencies results from the influence of inelastic 
scattering. For these low frequencies, the conductivity 
peak marking the onset of inelastic scattering appears 
just beyond the upper bound of our temperature range 
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FIG. 6. Fits to experiment of thermal regime microwave 
conductivity plotted (a) versus temperature for each of five 
frequencies and (b) versus frequency for each of four temper- 
atures. Points denote the temperature-dependent part of the 
microwave r«onductivity measured in YBa2Cu30e.993 by Hos - 
seini el alB Lines denote a two parameter fit of Eq. (5.5C) 
obtained with transport lifetime, Tt r 



2.93 x 10 _ii s, and 




Fermi-liquid renormalization factor, a fl = O.i 



FIG. 7. Results of exact numerical calculation of mi- 
crowave conductivity plotted as a function of temperature for 
(a) 1.14 GHz, (b) 13.4 GHz, (c) 22.7 GHz, and (d) 75.3 GHz. 
Straight lines connecting the first and last data points have 
been included to highlight any deviation from linearity. 



of interest. Quite naturally, the onset of inelastic scat- 
tering decreases the conductivity and yields a concave- 
down deviation from linearity. As we are concerned here 
with the nature of the low temperature elastic scatter- 
ing mechanism, we do not expect to reproduce this fea- 
ture. However, the observed high frequency (concave- 
up) deviations from linearity, as well as the temperature- 
independent shifts, do appear to be features of the elastic 
scattering w hich w e seek to understand. 

Since Eq. ( |5.5C| ) is an approximate result, obtained by 
taking the limit 1 /r, SI <C T and thereby neglecting sub- 
dominant terms, it is possible that sub-dominant features 
(constant shift and deviation from linearity for high fre- 
quencies) were neglected when we assumed the thermal 
limit to obtain our analytical result. To see if these fea- 
tures emerge in an exact solution, we present numerical 
results, valid beyond the thermal limit, in the following 
section. 



B. Numerical Results 



Recall from Sec. IV that the self-energy in the pres- 



ence of line defect scatter ing i s determined by the self- 
consistent solution of Eq. ( 4.13 ) as a function of momen- 
tum and energy. While an analytic expression was ob- 
tained in the thermal limit, this equation must be solved 
numerically for more general parameter values. Doing 
so (via Newton's method) yields the self-energy func- 
tions, S(p, m) a nd Si( p, uj). Pl uggin g these func tions 
into Eqs. (|5~29l), (|5.30|), (|]||), (0|), and (0|) and 



numerically integrating over momentum and energy, we 
obtain the microwave conductivity for a particular tem- 
perature and frequency. Repeating this process for all 
temperatures and frequencies of interest yields an exact 
(numerical) result for the microwave conductivity. 

Following this procedure, we computed the microwave 
conductivity as a function of temperature, from to 15 K, 
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for each of the five experimentally relevant frequencies. 
In all cases, we used values of the transport lifetime and 
Fermi-liquid correction factor obtained, in the previous 
section, from the fit of our thermal regime expression 
to experiment. The results for 1.14 GHz, 13.4 GHz, 
22.7 GHz, and 75.3 GHz are plotted versus temperature 
in Fig. 0. In each plot, a straight line has been drawn 
between the first and last data points to highlight any 
deviation from linearity. These exact results confirm the 
conclusions of our thermal regime calculation, revealing 
a predominantly linear temperature-dependence and a 
near-Drude frequency dependence. In addition, at the 
high frequencies where our thermal regime approxima- 
tions were least justified, we obtain, in qualitative agree- 
ment with experiment, concave-up deviations from lin- 
earity. Furthermore, these results do exhibit nonzero, 
albeit very small, offsets at zero temperature (most evi- 
dent for 75.3 GHz). However, the computed offsets are 
far smaller than those observed experimentally. Perhaps 
the point defect scattering, which is certainly present but 
was ignored herein for simplicity, is required to reproduce 
this feature correctly. A more complete calculation must 
consider the effects of both extended linear defects and 
point defects in the same system. This is an important 
direction for future study. 



VII. CONCLUSIONS 

The low temperature behavior of the microwave co 
ductivity measured in YBa2Cu3 06.993 by Hosseini et al 
has been difficult to explain in terms of point defect (im- 
purity) scattering. The basic problem is that the strong 
energy dependence of the low temperature quasiparticle 
density of states is directly reflected in a point defect 
scattering rate which is too energy dependent to produce 
a microwave conductivity that agrees with experiment. 
This suggests that an additional low temperature scat- 
tering mechanism may also be important. 

We note that extended linear defect scattering is an 
appealing candidate. Unlike point defect scattering, line 
defect scattering does not sample the full quasiparticle 
density of states, so the line defect scattering rate need 
not inherit a strong energy dependence. Furthermore, 
lines are prevalent in as-grown YBCO in the form of twin 
boundaries. Although these twin boundaries are elimi- 
nated via the application of a uniaxial stress, if remnants 
of the twinning structure are left behind, they could take 
the form of line defects. 

We have calculated (within the Born approximation) 
the self-energy and microwave conductivity due to scat- 
tering from a single domain of parallel, randomly spaced 
line defects, all aligned at —45° to the crystal axes. We 
find that the anisotropic nature of line defect scattering 
clearly differentiates odd nodes (for which electrical cur- 
rent is perpendicular to the defect lines) from even nodes 
(for which electrical current is parallel to defect lines). 



For the former, both intra-node (forward) scattering and 
opposite-node (back) scattering are permitted, whereas 
for the latter, only intra-node (forward) scattering is al- 
lowed. By including vertex corrections in our calcula- 
tion, we have accounted for the fact that back scattering 
is an effective means of degrading a current while for- 
ward scattering is not. Therefore, while odd nodes yield 
a nonzero transport scattering rate and microwave con- 
ductivity, even nodes make no contribution. The result- 
ing conductivity tensor is anisotropic, reflective of the 
anisotropic nature of the scattering mechanism. Thus, 
for a sample in which all twin boundaries were parallel 
prior to detwinning, this anisotropy should be observable. 
(In fact, the measurement of the conductivity tensor in 
such a single-domain sample would serve as a good test 
for the presence of line defects.) However, since the sam- 
ples in question contain multiple domains with line de- 
fects aligned at ±45°, we average over oppositely aligned 
domains to recover a conductivity scalar. 

In the experimentally relevant limit, 1/t, <C T <C 
Ao, which we have called the thermal regime, our calcu- 
lations simplify and we obtain an analytic expression for 
the microwave conductivity: 



a(fi,T) = 4-21n(2)/(f7r° r ) 

7T V 2 



1 



ikr ksT1 - {7A) 

As this result exhibits a linear temperature dependence 
and a near-Drude lineshape (modified loga rithmically 
by the function f(x) defined in Eq. ( 5.43| )), it cap- 
tures the most robust qualitative features of the Hosseini 
data. Furthermore, when we fit this expression to the 
temperature-dependent part of the measured conductiv- 
ity (see Sec. VIA), we obtain good quantitative agree- 
ment with reasonable values of our two fitting parame- 
ters, the transport lifetime and the Fermi-liquid renor- 
malization factor. However, we note that this approxi- 
mate analytical result fails to reproduce the more subtle 
features of the measured data: temperature-independent 
offsets and deviations of the temperature dependence 
from strict linearity at high frequencies. (We note that 
the deviations observed at low frequencies can be ex- 
plained by the onset of inelastic scattering.) 

In search of these subdominant features, we performed 
an exact numerical calculation, valid beyon d the thermal 
regime. These numerical results (see Sec. VI B) do ex- 
hibit deviations from linear temperature dependence at 
high frequencies, in qualitative agreement with experi- 
ment. While small temperature-independent offsets are 
also obtained, the magnitude of these is far smaller than 
observed experimentally. Perhaps the point defect scat- 
tering, which we have neglected herein, must be included 
to reproduce this feature. 

The close agreement of our calculated conductivity 
with experiment strongly suggests that extended linear 
defects are indeed present in detwinned single crystals of 
YBa2Cu3 6.993 and make an important contribution to 
the scattering. Our picture has been that these line de- 
fects are remnants of the twin boundary structure of the 
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as-grown crystal, left behind after detwinning. We de- 
scribed, in Sec. a detwinning scenario whereby the 
annihilation of twin boundaries leaves behind lines of 
oxygen vacancies in the CuO layer. While this is one 
example of a process that could yield extended linear de- 
fects, others are certainly possible. Our results, however, 
imply that something similar to this is taking place. 

To answer remaining questions and provide a more 
complete picture of the low temperature scattering, fu- 
ture calculations should include the effects of both line 
defect scattering and point defect scattering in the same 
system. Nevertheless, these initial results suggest that 
scattering from extended linear defects has a significant 
influence on the observed microwave conductivity in de- 
twinned single crystals of YBa2Cu3 06.993- 
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